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COSMOLOGICAL PARAMETERS FROM THE 2003 FLIGHT OF BOOMERANG 
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ABSTRACT 

We present the cosmological parameters from the CMB intensity and polarization power spectra of 
the 2003 Antarctic flight of the BOOMERANG telescope. The BOOMERANG data alone constrains 
the parameters of the ACDM model remarkably well and is consistent with constraints from a multi- 
experiment combined CMB data set. We add LSS data from the 2dF and SDSS redshift surveys to the 
combined CMB data set and test several extensions to the standard model including: running of the 
spectral index, curvature, tensor modes, the effect of massive neutrinos, and an effective equation of state 
for dark energy. We also include an analysis of constraints to a model which allows a CDM isocurvature 
admixture. 

Subject headings: cosmology, cosmic microwave background, polarization 



1. INTRODUCTION 

The angular power spectra of the Cosmic Microwave 
Background (CMB) have become invaluable observables 
for constraining cosmological models. The position and 
amplitude of the peaks and dips of the CMB spectra are 
sensitive to such parameters as the geometry of the Uni- 
verse, the cosmological constant and the energy densities 
associated with baryons and cold dark matter (Bond et al. 
1997). The CMB intensity spectrum has been measured 
with high precision on large angular scales (£ < 600) by the 
WMAP experiment (Hinshaw et al. 2003), while smaller 
angular scales have been probed by ground and balloon- 
based CMB experiments (Ruhl et al. 2003; Readhead et al. 
2004a; Dickinson et al. 2004; Kuo et al. 2004; Halverson 
et al. 2002). These data are broadly consistent with a 
ACDM model in which the Universe is spatially flat and 
is composed of radiation, baryons, neutrinos and the ever 
mysterious duo, cold dark matter and dark energy. 
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One of the firm predictions of this standard model is that 
the CMB is intrinsically polarized. Observations of the 
polarization power spectra, and the correlation with the 
total intensity spectra can therefore be used as a powerful 
consistency check, as well as potentially providing addi- 
tional cosmological information. On large angular scales 
the polarization is sensitive to the details of the reioniza- 
tion history and the curl component is a unique signature 
of tensor perturbations. On smaller angular scales the po- 
larization spectra can verify some of the basic assumptions 
made in the standard model. For instance, peaks in the 
polarization spectra arise from the same acoustic oscilla- 
tions at last scattering as those in the total intensity spec- 
tra. However, the peaks in the polarization spectra are 
predicted to be out of phase with the intensity peaks since 
the former are sourced by the velocity term of the photon- 
baryon fluid as opposed to its density. This effect provides 
the strongest constraint on the origin of the structure ob- 
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served in the spectra and breaks the severe degeneracy that 
is introduced in models with radically broken scale invari- 
ance. These are models in which non-trivial structure may 
already exist in the spectrum of initial perturbations. 

The recent polarization measurements of the DASI (Ko- 
vac et al. 2002), CAPMAP (Hedman ct al. 2002), WMAP 
(Kogut et al. 2003), and CBI (Readhead et al. 2004b) ex- 
periments have confirmed that the CMB is indeed polar- 
ized, providing an independent means for testing the un- 
derlying model. Many of the standard model cosmologi- 
cal parameters are becoming highly constrained, especially 
in combination with complementary data sets (e.g. Seljak 
et al. 2004). 

In this paper we test the standard model against the 
data from the 2003 long-duration balloon (LDB) flight of 
the BOOMERANG experiment (hereafter B03). This mis- 
sion marks the instrument's second successful trip over the 
Antarctic continent. The first LDB flight in December of 
1998 (hereafter B98), resulted in landmark, high signal- 
to-noise maps of the CMB intensity anisotropy and the 
detection of the first few peaks of the intensity angular 
power spectrum (de Bernardis et al. 2000; Netterfield et al. 
2002; Ruhl et al. 2003). For the 2003 flight the instrument 
receiver was redesigned to be polarization sensitive and 
the pointing system was upgraded to enable better atti- 
tude reconstruction. The B03 sky coverage is comprised 
of ^195 hours of data over <~1.8% of the sky, with an effec- 
tive beam 11.5 ± 0.23 arcminutes. Instrument calibration 
is based on the 90 GHz and 60 GHz WMAP data and the 
resulting amplitude uncertainty in calibration is ~ 2%. A 
complete instrument description, and the B03 CMB and 
Galactic maps are given in Masi et al. (2005). The final 
data set from the flight is comprised of four power spec- 
tra: the intensity power spectrum, TT; the EE (curl-free) 
and BB (curl-like) polarization power spectra; and the TE 
cross-power spectrum. These spectral data are presented 
in Jones et al. (2005b), Montroy et al. (2005) and Piacen- 
tini et al. (2005). 

This analysis examines in detail the cosmological impli- 
cations of the B03 data set. We begin by outlining our 
required data products and methodology in Section 2. We 
describe in Section 3 the various data combinations that 
are used in this analysis. In Section 4.1 we focus on the 
standard ACDM model and, applying only weakly restric- 
tive priors, we find that the simple parameter fits to B03 
data alone are fully consistent with those derived from 
other existing CMB data. To this CMB data, including 
the B03 data, we add in recent Large Scale Structure (LSS) 
rcdshift survey data, consisting of matter power spec- 
tra from the Sloan Digital Sky Survey (SDSS) (Tegmark 
et al. 2003) and the 2 Degree Field Galaxy Rcdshift Sur- 
vey (2dFGRS) (Percival et al. 2001), and determine the 
marginalized parameter constraints from this combined 
cosmological data set. In Section 4.2 we extend our analy- 
sis to include tests of several modifications of the standard 
model with the combined data sets. 

All of the models in Section 4 share the assumption that 
the initial perturbations of the primordial plasma are adi- 
abatic: in the early radiation dominated era the matter 
and radiation densities are all identically perturbed, giv- 

20 http://cmb.phys.cwru.edu/boomerang 
http: / /oberon. romal.infn.it/boomcrang/b2k 



ing an overall total density and hence curvature pertur- 
bation. This is not, however, the only possibility. Isocur- 
vature modes describe the other linear combinations of 
matter and radiation perturbations that do not contribute 
a curvature perturbation initially. Models with isocurva- 
ture contributions to the perturbations give rise to distinct 
signatures in the total intensity and polarization spectra. 
The latter can be used to further constrain the possible 
contributions by isocurvature modes that are not ruled out 
by measurements of the total intensity spectrum alone. In 
Section 5 we explore the constraints of the B03 and other 
data on a model with a mixture of a dominant adiabatic 
mode and a sub-dominant isocurvature mode. 

2. DATA PRODUCTS AND METHODOLOGY 

2.1. Summary of BOS Results 

We have developed two parallel and independent 
pipelines that we use to reduce the B03 observations 
from the time-ordered data to polarization and intensity 
anisotropy maps, through to angular power spectra. One 
pipeline was developed predominantly in North America 
(NA pipeline) (Contaldi et al. 2005; Jones et al. 2005a) 
and the other predominantly in Italy (IT pipeline) (Masi 
et al. 2005). The purpose of constructing two separate end- 
to-end pipelines is to check for self-consistency at various 
stages during the reduction and to check for robustness in 
the final spectra. We have carried out an extensive com- 
parison of the output of the NA and IT pipelines. We find 
excellent agreement for both the TT spectrum (Jones et al. 
2005b), representing the high signal-to-noise limit, and the 
polarization spectra (Montroy et al. 2005; Piacentini et al. 
2005) which are the most sensitive to the treatment of the 
experiment's noise characteristics and systematics. B03 
spectra obtained from the IT pipeline were tested with 
the same weak priors for the standard model case. The 
resulting parameter determinations are in good agreement 
with those reported here. 

The parameter constraints presented in this analysis are 
based on the output of the Xfaster hybrid Monte Carlo- 
maximum likelihood estimator (Contaldi et al. 20 05) 20 . 
The estimator uses a close to optimal, quadratic, Fisher 
matrix based estimator which is calibrated using signal 
only and noise only simulations of the entire data set, from 
time stream to final maps. It determines true polarization 
and total intensity angular power spectra (averaged over 
pre-determined £ bands) on the sky. After an arbitrary ini- 
tial guess, the quadratic estimator iterates onto the max- 
imum likelihood solution (Bond et al. 2000), C£ at , with 
errors determined by an estimate of the Fisher matrix for 
all band powers self-consistently. This ensures that the 
variance for each band power includes contributions from 
all cross terms and from all spectra. This is particularly 
important in the case of the cross spectra, as for example 
with the TE sample variance which is susceptible to the 
TT and EE power in addition to the TE power itself. 

The calculation of the full Fisher matrix also allows us 
to exclude band powers self consistently by cutting rows 
and columns from the inverse Fisher matrix. The effect of 
reduced sky coverage and/or pixel weighting is accounted 
for by computing all coupling kernels following Hivon ct al. 
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(2001) and Chon et al. (2004). The analysis typically in- 
cludes a simultaneous determination of a complete set of 
TT, EE, BB, TE, TB, and EB band powers. The EB 
and TB spectra are consistent with zero (as expected) and 
are excluded from the parameter determination by cutting 
out the bands in the inverse Fisher matrix (equivalent to 
marginalizing over their contribution). 

The spectra used in this analysis are shown in Figure 1. 
The data have been divided into bands which are generally 
A£ = 50 wide for TT and A£ = 100 wide for the three re- 
maining spectra. The multipole ranges for the B03 spectra 
which are used in this analysis are presented in Table 1. All 
band-to-band correlations are included in the Fisher infor- 
mation matrix and are at most ^20%. The band spacing 
was chosen in part to ensure that these correlations were 
not large. 

The Xfaster code also calculates the required band 
window functions, Wf, which are used to convert the 
model power spectra, C™ od , into theoretical bandpowers 
via 

i[wfcr d ] (1) 



<cr d > 



mi 



Here Cf od 



£(£ + l)C7 od /27r and we have introduced 
the notation for the "logarithmic integral" of a spectrum 

(Bond et al. 2000), l[f t ] = J2e m+tjf^ Thc above opera- 
tion permits direct comparison of theory C B od with data 

B ' 

A final issue is the potential bias introduced by the non- 
Gaussian distribution of thc bandpowers in the signal dom- 
inated regime. It has been shown (Bond et al. 2000) that 
the variable Z B = ln(C^ at + Cg ) is more normally dis- 
tributed than thc bandpowers C B . The noise offsets, C B , 
are a measure of the deconvolved noise spectrum on the 
sky and are calculated with the same quadratic estima- 
tor using Xfaster on the average of simulated noise-only 
observations. 

Thc distribution of the bandpowers tends to a Gaussian 
in the noise dominated regime and log-normal in thc sam- 
ple variance dominated regime. Both limits are significant 
for thc TT bandpowers, hence we transform all the TT 
bands to offset log-normal variables and treat the likeli- 
hood function in the new variables as Gaussian for param- 
eter estimation. For the polarization spectra EE and BB, 
which are noise-dominated, we use Zb = C B , with no 
non-Gaussian correction. For TE we also use Zb — C B at 
since negative values of C B at occur. The Fisher matrix 
of the bandpowers is transformed as Fbb' = Z' b Fbb'Z' b , 
with Z' B = dZ B /dC B at = (Cf 1 * + Cg)- 1 if B is a TT 
bandpower and Z' B = 1 otherwise. 

In summary the Xfaster data products include the 
bandpowers, Fisher matrix, window functions and noise 
offsets. 

2.2. Parameter Estimation Methodology 

Thc Monte Carlo Markov Chain (MCMC) sampling 
technique we use for parameter estimation is described 
in detail in Neal (1993); Christensen et al. (2001); Lewis 
& Bridle (2002) and implemented in the publicly avail- 
able COSMOMC 21 package. Here we give a brief summary 

21 http://cosmologist.info/cosmomc 



of the relevant details. Thc technique uses a Bayesian 
approach, generating samples of the posterior probability 
density function (PDF) of the parameters y given the data 

P(y\z) ex P(y)P(z\y), (2) 

where P(z\y) is the likelihood PDF and P(y) is the prior 
PDF of y. The posterior is sampled by running a num- 
ber of Markov Chains. The chains are constructed via the 
Metropolis-Hastings (MH) algorithm whereby a candidate 
parameter vector, y' is determined from an arbitrary pro- 
posal density distribution q(y'\y n ) where y n is the current 
state of the chain. The candidate y' is accepted with ac- 
ceptance probability given by 



• , P(y'\z)g(y n \y') 



(3) 



At each point in the chain the acceptance probability for a 
candidate point is compared to a random number u drawn 
uniformly in thc to 1 range. If u < ct(y'\y n ) then the 
proposed vector is accepted and the next point in the chain 
is y n +i = y' ■ If u > a{y'\y n ) then the proposed vector 
is rejected and y n +i = y n - 

For the B03 CMB data the likelihood evaluation at each 
point in the chain requires the calculation of 

x 2 = E ( z B od (y) - z d B at )F B B> (z^ d (y) - zp). (4) 

BB' 

The WMAP data likelihood is computed using the like- 
lihood code supplied by the WMAP team (Verde et al. 
2003; Kogut et al. 2003), but with two modifications. The 
first modification is a change to the TE likelihood func- 
tion to account for the correlation between the intensity 
and TE power spectrum estimators (we neglect the small 
correlations between the Ci estimators at different I) (Dore 
et al. 2004). After the chains have been run we use impor- 
tance sampling (e.g. see Lewis & Bridle (2002)) to correct 
the WMAP likelihood on large scales using the more com- 
putationally intensive likelihood code from Slosar et al. 
(2004). This Slosar-Seljak modification uses a more accu- 
rate calculation of the WMAP likelihood at low multipoles 
(£ < 11) and considers in more detail the errors associated 
with foreground removal. 

The theoretical CMB spectra (as well as the matter 
power spectra) are computed using CAMB (Lewis et al. 
2000), a fast parallel Boltzmann code based on CMB- 
FAST (Seljak & Zaldarriaga 1996). We calculate statistics 
of interest, such as the marginalized posterior distribution 
of individual parameters, from the MCMC samples after 
removing burn in. We run six chains for each combination 
of data and parameters that cannot be importance sam- 
pled. We marginalize numerically over each data point's 
calibration and beam uncertainties at each sample in the 
chain. The calibration errors are assumed to be indepen- 
dent between data sets. We check convergence by ensuring 
that the standard deviation between chains of the 95%- 
pcrcentilc estimated from each chain is less than 0.2 in 
units of the all-chain parameter standard deviation. This 
should ensure that sampling errors on quoted limits are 
minimal. 

Parameter estimates from MCMC have been shown to 
be in very good agreement with those derived using an 
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Fig. 1. — The B03 bandpowers used in this analysis. We have included the total intensity TT, polarization EE and BB, and cross 
correlation TE spectra. The EB and TB spectra are excluded from this parameter analysis. The solid/black curve is the previous 
concordance model, a best fit to WMAP (first-year) +CBI+ACBAR data from http://lambda.gsfc.nasa.gov/product/map/, with 
(Q, b h 2 ,n c }i 2 ,n s (k = 0.05),exp(-2r), A(fc = 0.05), h) = (0.0224, 0.111, 0.958, 0.802, 0.739, 0.720). The yellow/dotted curve is the 
CMBall (Table 3)+B03 maximum likelihood ACDM model from this analysis with (slightly different parameterization-see text), 
(n b h 2 ,n c h 2 ,n s (k = 0.05), -r,ln(10 10 A s (fc = 0.05)), 6) = (0.0228, 0.108, 0.959, 0.138, 3.12, 1.04). 



adaptive C^-grid (Bond et al. 2003) that was previously 
applied to the B98 analysis (Ruhl et al. 2003). MCMC 
results for CMBall+B98 (Bond et al. 2004) are also in 
good agreement with those we obtain for CMBall+B03 
for the baseline model defined below with the same priors 
applied 22 . 

3. DATA COMBINATIONS 

3.1. The CMB data 

We consider a number of combinations of data. We 
break the B03 data set into one subset consisting of the 
TT spectrum alone (B03TT), and another subset consist- 
ing of the EE, BB and TE spectra (B03pol) alone. We 
also consider fits to the entire B03 data set, WMAP data 
alone and a combined B03 + WMAP data set. We next 
combine B03 with available data from a collection of CMB 
experiments. We outline in Table 3 the experiments and 
multipole ranges which make up that collection, which we 

22 B03 and B98, with overlapping sky coverage, are correlated data 
combined B98 and B03 maps in a future analysis. 



call CMBall. We note that because of the overlap in I 
range of the ARCHEOPS (Tristram et al. 2005) data with 
the WMAP data, the former cannot be included in the 
CBMall data set, unless a joint analysis is done. The B03 
multipole range is given in Table 1. The cosmic variance 
of the WMAP and B03 data sets is correlated in the low 
multipole range (essentially over the first peak of the TT 
power spectrum). To account for this, we cut the lower 
multipoles of the B03 TT spectrum (£ < 375) when com- 
bining B03 data with WMAP data. 

3.2. The LSS Data 

For our final data combination we also include LSS ob- 
servations from 2dFGRS and the SDSS. The two redshift 
surveys are treated in a conservative fashion. For exam- 
ple, although the SDSS bandpowers have been corrected 
for differing galaxy bias factors associated with different 
types of galaxies, there is still an overall galaxy bias fac- 
tor, b g , the ratio of the square root of the galaxy-galaxy 

We therefore exclude B98 from this analysis and will consider the 
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power spectrum for galaxies to that of the mass den- 
sity power spectrum today. Although the indications are 
that this is a number near unity (Percival et al. 2001; 
Tegmark ct al. 2003), in our standard results we allow 
it to take on arbitrary values by marginalizing it over a 
very broad distribution. This means that our LSS infor- 
mation is only constraining models through the shape of 
the power spectrum, but not the overall amplitude. Con- 
straining the overall amplitude is akin to imposing a prior 
on as- To test sensitivity to this, we have adopted var- 
ied Gaussian errors on b 2 about a mean. We have taken 
the mean to be unity and adopted errors on Sb 2 appro- 
priate for 5b g = 0,0.1,0.5 and 10, then marginalized over 
b 2 g . A uniform prior in b 2 g leads to the same results as for 
Sb g = 10 23 . Most parameter averages we obtain are rela- 
tively insensitive to Sb g . We comment on its effect below: 
it has impact on the massive neutrino and dark energy 
equation of state constraints. We only use SDSS data for 
wavenumbers k < O.lftMpc -1 to avoid nonlinear correc- 
tions and to avoid possible non-uniform b g complications. 
(See Tegmark et al. (2003) for a discussion of these and 
other issues.) A similar Sb g marginalization strategy was 
used for the 2dFGRS rcdshift survey data. 

An estimate using galaxy-galaxy lensing from SDSS 
(Seljak et al. 2005) is b g = 0.99 ± 0.07. (These authors 
also used WMAP data to obtain this value, so it is not a 
completely independent determination of the bias.) An 
estimate using the 3-point function and rcdshift space 
clustering distortions for 2dFGRS gives b g = 1.04 ± 0.04 
(Verde et al. 2002). Based on these two analyses, adopting 
b g = 1.0 ± 0.10 to illustrate the effect of knowing the bias 
better, which translates into a as prior, seems reasonable. 

For the purposes of this paper, in which our focus is on 
the B03 CMB data, we have limited the LSS information 
we include. For example, we have not incorporated the 
SDSS results on luminous red galaxies (Eisenstein et al. 
2005). The recent final power spectrum and window func- 
tions of the 2dFGRS survey Cole et al. (2005) is not yet 
available. 

3.3. Other Data Sets 

We have applied the supernova data (SNIa) in Sec- 
tion 4.2.5 to the determine the dark energy equation of 
state. For this we use the gold set, as described in Riess 
et al. (2004). Also, for a few cases we include the Hq prior 
value from the HST Key Project (Freedman et al. 2001). 

We do not explicitly include weak lensing results. These 
generally determine the parameter combination as^l^ 
providing additional independent constraints (Contaldi 
et al. 2003; Bond et al. 2005; Readhead et al. 2004b). We 
also do not include information on the Lyman alpha for- 
est, even though it probes the power spectrum to smaller 
scales. Although adding this data does result in some more 
stringent constraints than those we derive here (McDonald 
et al. 2004), the forest information is more susceptible to 
scale dependent biasing effects associated with gasdynam- 
ical and radiation processes. 

4. ADIABATIC MODELS 

4.1. Baseline Model 

23 Note that allowing b 2 to be negative has no effect and yields the same 



4.1.1. Parameterization and Priors 

For our baseline model we consider a flat universe with 
photons, baryons, massless neutrinos, cold dark matter 
and a cosmological constant. Initial conditions will be 
taken to be purely adiabatic (no isocurvature modes). We 
assume a power law form for the power spectrum of the 
primordial comoving curvature perturbation, described by 
V s = A s (k/K)( n °-V, where the n s is the scalar spectral 
spectral index and A s is the scalar amplitude (we choose 
a pivot point k* = 0.05Mpc _1 ). The physical baryon den- 
sity and dark matter density are parameterized by fli,h 2 
and fl c h 2 , where h = i/ /100km s _1 Mpc _1 is the Hub- 
ble parameter. We use the parameter to character- 
ize the positions of the peaks in the angular power spec- 
tra, defined as one hundred times the ratio of the sound 
horizon to the angular diameter distance at last scatter- 
ing (Kosowsky et al. 2002). Finally, the parameter r is 
used to describe the Thomson scattering optical depth to 
decoupling. Thus our baseline model is a function of 6 
cosmological parameters to which we impose the following 
flat weak priors: 0.5 < n s < 1.5; 2.7 < ln(10 10 A s ) < 4.0; 
0.005 < n b h 2 < 0.1; 0.01 < n c h 2 < 0.99; 0.5 < 6 < 10.0; 
and 0.01 < r < 0.8. Additional weak priors restrict the 
age of the universe to lOGyr < age < 20Gyr and the 
expansion rate to 0.45 < h < 0.9. All priors are sum- 
marized in Table 2. Besides being generally agreed upon 
by cosmologists, our weak priors are consistent with those 
used in much of the CMB literature, e.g., Lange et al. 
(2001), Bond et al. (2003) and Readhead et al. (2004b). 
We choose not to impose the restrictive prior t < 0.3 ap- 
plied in Spcrgcl et al. (2003). We note that some of our 
results are sensitive to our choice of prior on H and we 
explore the effect of strengthening our Hq prior in Sections 
4.2.1 and 4.2.2. 

In addition to the base parameter values, the results 
also include marginalized constraints for several derived 
parameters including: Q\, the relative dark energy den- 
sity; the age of the Universe; O m , the relative total matter 
density; as, the root mean square linear mass perturba- 
tion in 8/i _1 Mpc spheres; z re , the redshift of reionization 
assuming it is a sharp transition; and the Hubble constant, 

4.1.2. Consistency of BOS Data Set 

The resulting marginalized parameter constraints for 
the baseline model for each of the data combinations are 
given in Table 4 and presented graphically in Figure 2. In 
both Table 4 (and in the ones that follow) and Figure 2 
we give the Bayesian 50% probability value (the median) 
obtained from the marginalized probability for each pa- 
rameter. The quoted errors represent the 68% confidence 
interval obtained by integrating the marginalized distri- 
butions. In the case of upper or lower bounds, the 95% 
confidence limits are quoted. We note that our baseline 
CMBall+B03+LSS result is fairly insensitive to Sb g and 
that we have chosen the less restrictive flat, uniform prior 
in b 2 g . 

The comparison of B03pol and B03TT provides a robust 
internal consistency check. We note that the B03pol con- 
straints to flbh 2 and fl c h 2 are quite good with uncertain- 
ties which are only slightly larger than those of the B03TT 

results as a (uniform) positive b 2 constraint. 
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Fig. 2. — Median values obtained from the marginalized probability for each parameter for the baseline, standard model. The errors bars 
represent the 68% confidence interval. The 95% upper limit is given for the case of r for B03 data alone. The following flat weak priors are 
imposed (as outlined in Table 2): 0.5 < n s < 1.5; 2.7 < ln(10 & A s ) < 4.0; 0.005 < Q b h 2 < 0.1; 0.01 < Q c h 2 < 0.99; 0.5 < 6 < 10.0; and 
0.01 < t < 0.8. Additional weak priors restrict the age of the universe to lOGyr < age < 20Gyr and the expansion rate to 45 < Ho < 90. 
Our baseline CMBall+B03+LSS result is fairly insensitive to Sb g and we have chosen for this case the less restrictive flat, uniform prior in b 2 . 



result. However, the B03pol constraints on n s , r and A s combined parameter A s e~ 2T . The latter determines the 
are weak and results for these cases are prior driven. We overall power in the observed CMB anisotropy (except at 
present in Figure 3 a 2D likelihood plot of 9 versus the low £), and is therefore better constrained than the pri- 
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mordial power ^4 S . CosmoMC uses a covariance matrix 
for the parameters and is therefore able to ascertain lin- 
ear combination degeneracies. Although we use In A, and 
t as base parameters, the proposal density knows that 
the combination ln(A s e T ) is well constrained and can 
explore the poorly constrained orthogonal direction effi- 
ciently. We find that the B03 data alone does particularly 
well at constraining A s e~ 2r . The angular-diameter dis- 
tance variable 8 defines the shift with £ of the overall Cg 
pattern, in particular of the pattern of peaks and troughs. 
With all of the CMB data it is the best determined pa- 
rameter in cosmology, 1.045 ± 0.004; with B03pol it is an 
important test which demonstrates the consistency of the 
positions of the polarization spectra peaks and troughs 
relative to those forecasted from the TT data, although 
the errors are larger with 6 = 1.08 ± 0.03. For the CBI 
TT, TE and EE data in combination with WMAP TT 
and TE, Readhead et al. (2004b) found 6 = 1.044 ±0.005. 
With just the CBI EE polarization data they determined 
6 = 1.06 ± 0.04, again showing the consistency we find of 
the data with the TT forecast of the polarization peaks 
and trough. 

The B03 median parameter values are remarkably con- 
sistent with the parameter constraints from WMAP data 
alone. We note that in general the Slosar-Seljak modi- 
fication to WMAP tends to broaden WMAP parameter 
likelihood curves and that the most significant impact on 
the median values is in r (<~ 0.3cr increase) and in Q m 
(~ 0.6cr decrease). Adding B03 to the WMAP data de- 
creases the parameter uncertainties by an average of 16%. 
The most significant effect is a ^33% decrease in the as 
uncertainty. Figure 4 shows the likelihood curves for the 
6 base parameters and 6 derived parameters for a variety 
of data combinations. Overall the various data combina- 
tions are generally in good agreement at better than the 
la level. The largest outlier is fl c h 2 which increases by 
1.5(7 with the addition of the LSS data set. Also, similar 
to Spcrgcl et al. (2003), we find that the addition of small 
scale CMB data lowers both the value for the amplitude of 
fluctuations at k — O.OSMpc^ 1 and the value of the scalar 
spectral index. The effect of adding the LSS data follows 
this trend. 

4.2. Modified Standard Model 

In this section we explore five extensions of the standard 
model by adding, in turn, one parameter to the baseline 
parameter set. In all cases we maintain the same weak pri- 
ors on the base parameters as outlined in Table 2. Some 
of the results are sensitive to our chosen prior range for 
H . For example, in certain cases we will note the impact 
of strengthening our H prior to the value from the HST 
Key Project (Freedman et al. 2001), h = 0.72 ±0.08, with 
the errors treated as Gaussian. 

4.2.1. Running Index 

We modify the power law form for the power spectrum 
of the density perturbations to allow the spectral index, n Sl 
to vary with scale. Following Kosowsky & Turner (1995) 
this variation can be parameterized by the term n run = 
dn s /dlnk, such that n s = n s (k+) + n lun (k+) ln(k/k+), 
where again fc* = 0.05Mpc _1 . We restrict n run to lie 
between -0.3 and 0.3. Results from the combined data 



sets, CMBall+B03 and CMBall±B03±LSS, are given in 
Table 5. 

Spcrgcl et al. (2003) report a detection of the run- 
ning index of n run = — 0.031^q'q]^ from their combined 
WMAPext+2dFGRS+Lyman a data set. Slosar et al. 
(2004) present a reduction in significance of the detection 
of 7i run when their full likelihood analysis and detailed fore- 
ground removal is applied to the WMAP data. We find 
that the Slosar-Seljak modification to WMAP decreases 
the significance of n run , but that inclusion of the data 
from the small scale CMB experiments has the opposite 
effect (as was the case found by Spergel et al. (2003)). 
From CMB data alone we determine a median value for 
n run = —0.072 ± 0.036. This result is somewhat sensi- 
tive to our choice of prior. Spergel et al. (2003) apply 
a strong t < 0.3 prior which effectively reduces the me- 
dian value of n run for their CMB data only case. Here 
we apply a Gaussian HST prior on Hq which lowers the 
significance of the running index to n run = — 0. 065^034 
for the CMBall + B03 data set. Inclusion of the LSS data 
(with uniform prior in b 2 ) further reduces the significance 
and our median value from the larger combined data set 
is n run = — 0.05lto;o26- We note that the application of 
a Gaussian prior to b 2 has no impact on the running in- 
dex parameter. Application of the HST prior on H yields 
a final median value n lun = —0.048 ± 0.026 for the CM- 
Ball + B03 + LSS (+HST) data set. Figure 5 shows the 
likelihood curves for the n run parameter for various data 
combinations. It is interesting to compare our result with 
that of Seljak et al. (2004), who argue that if the state-of- 
the-art modeling of Lyman a forest measurements is dom- 
inated by statistical rather than systematic errors, then 
|«run| < 0.01 

4.2.2. Curvature 

We consider a modification to the standard model which 
allows the possibility of non-flat geometry. We parame- 
terize the curvature density by Clk and allow it to vary 
between -0.3 to 0.3. Table 5 shows the results for the CM- 
Ball+B03 and CMBall+B03+LSS data sets. The CMB 
data alone places a constraint on the curvature which is 
n k = -0.037^-^. We show in Figure 6 the likelihood 
profiles for WMAP, WMAP±B03, CMBall±B03 and CM- 
Ball+B03+LSS. While the addition of B03 data to the 
WMAP data tends to lower the significance of curvature, 
adding more small scale CMB data increases the width of 
the low end tail. Addition of the LSS data, with uniform 
prior in b 2 , yields a median value of fifc = —0.027 ± 0.016. 
Application of the Gaussian prior in b 2 (with 10% uncer- 
tainty in b g ) has a slight effect with a resulting median 
value of Slfc = —0.022 ± 0.015. If we restrict the Hq value 
by the application of a Gaussian HST prior, the curvature 
density determined from the CMBall + B03 data set is 
Slfc = —0.015 ± 0.016. Moreover, application of the more 
stringent Hq prior reduces the median value of the curva- 
ture from the combined CMBall + B03 + LSS data set 
(flat b 2 prior) to Q k = -0.025 ± 0.017. Our result agrees 
well with the constraint = —0.010 ± 0.009 obtained by 
combining CMB data with the red luminous galaxy clus- 
tering data, which has its own signature of baryon acoustic 
oscillations (Eisenstein et al. 2005). 
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Fig. 3. — Constraints on A B e~ 2T versus 6. Inner contours represent 68% likelihood regions and outer contours 95% likelihood 
regions. The peak position characterization parameter 9 is best the determined parameter in cosmology, 1.045 ± 0.004, from 
the CMBall+B03 data set. We find that the B03TT data does particularly well at constraining both the peak pattern and the 
combined A B e~ 2T amplitude parameter. The constraint from WMAP alone on A B is better than that from B03. The agreement 
between the B03pol and B03TT data is consistent with the basic inflation picture. 



4.2.3. Tensor Modes 

So far we have assumed only scalar perturbations. How- 
ever inflationary models can produce tensor perturbations 
from gravitational waves that are predicted to evolve in- 
dependently of the scalar perturbations, with an uncor- 
rected power spectrum Vt- The amplitude of a tensor 
mode falls off rapidly after horizon crossing and the effect 
is therefore predominantly on the largest scales: tensor 
modes entering the horizon along the line of site to last 
scattering distort the photon propagation and generate an 
additional anisotropy pattern. We parameterize the tensor 
component by the ratio A t /A s , where A t is the primordial 
power in the transverse traceless part of the metric tensor 
on O.OSMpc^ 1 scales. We impose a very weak prior on the 
amplitude ratio, restricting it to lie between to 20. 

A tensor spectral index, defined by Vt oc k nt , must also 
be set. In inflation models it is related to the amplitude 
ratio by A t /A s as — 8n t /(l — n t /2), so one parameter suf- 
fices. In a nearly uniformly accelerating regime n t ~ n s — 1 
is also expected. However, although n s > 1 can arise in 
inflation models, n t > is difficult to obtain. As well, 
for many inflation models with n s — 1 just below zero, n% 
comes out to be slightly closer to zero. If the acceleration 
changes significantly over the observably range, n t and n s 
would not be intimately tied. Rather than let nt float as a 
second added parameter, we have chosen to make Vt flat 



in k (and thus set n t to zero) for the computations of the 
tensor-induced component of C(. 

Results are presented in Table 5, and Figure 7 illustrates 
the likelihood curves for the amplitude ratio for a number 
of data combinations. The influence of the high precision 
of the WMAP data on the largest scales is evident. Adding 
the small scale CMB data only slightly reduces the limit. 
We determine an upper limit on the tensor ratio from CMB 
data (CMBall+B03 data set) alone of A t /A s < 0.71 (95% 
confidence limit). The CMB data appear to select mod- 
els with relatively large tensor-to-scalar ratios. However, 
these models which have large values for the Hubble pa- 
rameter (Hq ~ 85) are allowed due to the poor constraint 
on Hq when including tensor modes. In this case, the con- 
straints on Hq are driven mainly by our choice of weak 
priors and the data only provides a lower limit (see Ta- 
ble 5) . We find with application of the HST prior (which 
excludes these models with large Hq values) that the ten- 
sor limit from the CMBall+B03 data set is reduced to 
At /A s < 0.635. A similar effect is obtained with the ad- 
dition of the LSS data which further reduces the limit to 
At/A s < 0.36. When we constrain b g in the LSS data, 
the limits are very similar. The application of the more 
restrictive prior discussed above, with only n s < 1 allowed 
to have a tensor contribution, lowers the CMBall+B03 
limit to A t /A s < 0.45 and the CMBall+B03+LSS limit 
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Fig. 4. — Marginalized one-dimensional distributions for the baseline model parameters for the data combinations WMAP only 
(black/dotted), WMAP + B03 (green/solid), CMBall + B03 (blue/dashed), and CMBall + B03 + LSS (red/dash-dotted). The curves 
are each normalized by their peak values. All distributions are derived from chains run with the weak set of external, uniform priors shown 
in Table 2. The LSS data consists of the 2dFGR,S and SDSS redshift surveys (with a flat bg prior imposed). The most significant impact of 
the B03 data is on erg. Moreover, the ag constraint from CMB data alone is quite strong, with the addition of LSS data having little effect. 



to A t /A s < 0.31. As a final case we set n t — —(At/A s )/8 
and found A t /A s < 0.54 for the CMBall+B03 data set and 
A t /A s < 0.30 for the combined CMBall+B03+LSS results 
(again applying the b g constraint has little effect). 

4.2.4. Massive Neutrinos 

Observational evidence from solar and atmospheric neu- 
trino experiments, such as the Sudbury Neutrino Observa- 
tory (Ahmad et al. 2002) and Super-Kamiokande (Toshito 
et al. 2001), suggest that neutrinos change flavour: neu- 
trinos of different generations oscillate into each other. 
The implication of flavour changing is that neutrinos have 
mass. Given that neutrinos are the second most abun- 
dant particles in the Universe, massive neutrinos could 
have considerable impact on the energy density of the 
early Universe. We consider here the case of three neutri- 



nos of degenerate mass, such that il„h 2 = 3m„/94.0eV. 
This assumption is well justified given the small square 
mass difference measured by oscillation experiments (at 
most 5ml ~ 10~ 3 eV, Aliani et al. (2003)). We param- 
eterize the massive neutrino contribution as a fraction of 
the dark matter energy density, /„ = Q u h 2 /SlDhih 2 = 

1 — ilcDAlh 2 /£lDNlh 2 . 

Results for the combined data sets are given in Table 5. 
From CMB data alone the upper limit on the neutrino frac- 
tion is /„ < 0.21 (95% confidence limit). This translates 
to an upper limit on the neutrino mass of m v < 1.0 eV 
or £l„h 2 < 0.033. This limit is more stringent than the 
3 eV upper limit on the electron neutrino mass determined 
from tritium beta decay experiments and recommended 
in the Review of Particle Physics (Eidelman et al. 2004). 
Including the LSS data (flat b 2 prior) pushes this limit 
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Fig. 5. — Marginalized one-dimensional distributions for the 
n IU n parameter for the baseline + running index model. Weak 
priors imposed are those outlined in Table 2. The running index 
parameter is restricted to lie between -0.3 and 0.3. Application 
of the HST prior on Ho slightly reduces the significance of a 
running index. 



Fig. 6. — Marginalized one-dimensional distributions for Qk 
for the baseline model which allows non-zero curvature. Weak 
priors imposed are those outlined in Table 2. We restrict f2fc 
to the range -0.3 and 0.3. The relatively wide scope for posi- 
tive curvature is associated with the angular-diameter-distance 
degeneracy which is only partly broken by the CMB data. Ap- 
plication of the HST prior on Ho to the larger combined data 
set somewhat reduces the possibility of significant curvature. 




Fig. 7. — Marginalized one-dimensional distributions for the 
amplitude ratio A t /A s for the baseline model modification 
which allows tensor modes. Weak priors imposed are those 
outlined in Table 2. We impose the weak prior < A t /A s < 20 
to the tensor contribution. We find from CMB data alone (CM- 
Ball+B03) an upper limit (95% confidence) on the amplitude 
ratio of A t /A s < 0.71. For these models however, Ho is only 
poorly constrained (see text). Addition of the LSS data reduces 
this limit to A t /A s < 0.36. 



Fig. 8. — Marginalized one-dimensional distributions for 
m v for the baseline model which allows massive neutrinos (3 
species of degenerate mass). Weak priors imposed are those 
outlined in Table 2. We parameterize the massive neutrino 
contribution as a fraction of the dark matter energy density, 
/„ = i} v h 2 /n DM h 2 . We find from CMB data alone (CM- 
Ball+B03) an upper limit (95% confidence) on the neutrino 
mass of m„ < 1.0 eV. Adding the LSS data reduces this limit to 
m„ < 0.40 eV, without any b g constraint, and to m v < 0.16 eV, 
when b g — 1.0 ± 0.10 is used. 
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down considerably to /„ < 0.093 (95% confidence) and 
limits the neutrino mass to m v < 0.40 eV. This result is 
somewhat larger than that found in Spergel et al. (2003). 
We find that addition of more and more small scale CMB 
data drives the limit up as is evident in Figure 8. When 
b g = 1.0 ± 0.10 is used, the neutrino fraction upper limit 
is reduced to /„ < 0.041 (95% confidence), correspond- 
ing to a neutrino mass limit of m v < 0.16 eV. This neu- 
trino mass limit is in good agreement with the strong limit 
(m v < 0.18 eV) obtained by (Seljak et al. 2005), who in- 
cluded the bias constraint and the SDSS and WMAP data. 
In their analysis of b g they found er 8 = 0.85+o;og, with 

b g = 1.02^0 08- This compares with the values we obtain: 
cr 8 = 0.85 ± 0.04 with Sb g = 0.10 and cr 8 = 0.74 ±0.08 with 
5b g = oo. 

4.2.5. Dark Energy 

The standard model predicts (and CMB observations 
strongly support) a universe which is nearly flat, implying 
a total energy density approaching critical. The total mat- 
ter density however, comprises only one third of the total 
energy density. The prevailing energy density component 
comes from some form of dark energy which up to now we 
have assumed takes the form of a vacuum density or cos- 
mological constant, A, with equation of state described by 
w = p/ p = — 1, where p and p arc the dark energy pressure 
and density respectively. We now consider the possibility 
that the dark energy component is a rolling scalar field or 
quintessence (see for example Ratra & Peebles (1988) or 
Hucy et al. (1999)), allowing the effective constant equa- 
tion of state parameter w to differ from — 1. We treat w as 
a redshift-independent phenomenological factor and allow 
it to range with a uniform prior over the range —4 to 0. We 
have also run the cases with w restricted to lie in the range 
— 1 to and find similar limits. To be self-consistent, per- 
turbations in the dark energy should be allowed for when 
to is not — 1 , although these have a small impact and only 
at low multipoles. We set the effective sound speed for 
the perturbations to unity in CAMB, the value for a scalar 
field. 

The marginalized one-dimensional distributions for var- 
ious data combinations are presented in Figure 9. We 
find from CMB data alone to = — 0.86±g;||. The addi- 
tion of the LSS data, applying the conservative uniform 



the degeneracy. 



5. SUB-DOMINANT ISOCURVATURE MODEL 

CMB anisotropics provide a powerful probe of the na- 
ture of early universe perturbations. However almost any 
TT power spectrum shape can be fit rather well by us- 
ing contrived combinations of initial perturbations, for ex- 
ample by adding structure to the primordial power spec- 
trum, and/or by adding isocurvature modes (which come 
in many varieties). With full freedom, determination of 
the basic cosmic parameters suffers because of high corre- 
lation with these extra degrees of freedom. These degen- 
eracies are broken by the addition of polarization data, be- 
cause of the pattern differences of the peaks and troughs 
between EE, TE and TT power spectra. In particular 
the intensity and polarization power spectra for isocurva- 
ture modes have the peaks out of phase with those from 
adiabatic modes (see e.g., Bond & Efstathiou (1987), Hu 
et al. (1997)). A mix of isocurvature and adiabatic modes 
can be designed to give acceptable fits to the CMB inten- 
sity power spectra (e.g., (Bucher et al. 2004; Kurki-Suonio 
et al. 2005)), and the polarization data, from B03 as well 
from DASI, CAPMAP and CBI are not yet at the point 
to clearly distinguish among these more complex models. 

To illustrate the constraints that can be determined 
from the current CMB data, we consider here a simple 
hybrid case consisting of our basic adiabatic mode model 
with constant spectral index, a single cold dark matter 
(CDM) isocurvature mode with its own constant primor- 
dial spectral index n; so , with no correlation between the 
two. This adds another two parameters to our basic six, 
n- lso and an amplitude ratio i?2 = {Ai ao /A s ). We assume 
the isocurvature perturbations are Gaussian-distributed as 
we have done for the adiabatic modes. Results are shown 
in Table 6 for the CMBall+B03±HST data combination. 
Aside from the more stringent HST data prior on Hq, all 
priors on the 6 base parameters are as outlined in Table 2. 
Although results indicate that there is no evidence for the 
presence of an isocurvature mode, the upper limits still 
allow for a sub-dominant component. 

We now expand on the theoretical framework. Isocur- 
vature modes may arise in two (or more) field models of 
inflation, and, though certainly not a natural prediction, 
they have reasonably good physical motivation. Isocurva- 
ture modes could also be generated after inflation ended. 
An ingredient needed for isocurvature modes to have an 
observable impact on the CMB is that they are associated 
with a component of significant mass-energy. If the dark 
matter is cold and of one type there are two distinct matter 
isocurvature modes, the baryon and CDM modes, involv- 
ing primordial fluctuations in the entropy-per-baryon or 
the entropy-per-CDM-particle. The classic example of a 
CDM possibility is the isocurvature axion mode. It turns 
out that the isocurvature CDM and baryon modes actually 
have almost identical signatures in the CMB (Gordon & 
Lewis 2003), and hence cannot be constrained separately. 
In this paper we choose to constrain only the CDM isocur- 
vature mode. There are other possibilities for isocurvature 
modes that we will not explore here 24 . 

24 One set involves neutrinos and photons compensating each other (Bucher et al. 2000; Rebhan & Schwarz 1994; Lewis 2004), but these must 
be produced after neutrino decoupling and so seem quite unlikely. Much better motivated are isocurvature modes associated with defects such 



flat prior on b 2 to the galaxy bias factor, yields a me- 
dian value of to = —0.65 ±0.16. This result is highly 
sensitive to our choice of prior on b 2 g . The uniform flat 
prior on b 2 g gives a relatively high best fit bias value of 



1.3. Applying a more restrictive Gaussian prior to 



a 
w 



-0.95 ± 0.15 with b„ 



1.0 ± 0.10; and 



gives: w 

= -0.75 ±0.15 with b g = 1.0 ±0.50. We explore the ef- 
fect of adding the SNIa data which significantly improves 
the constraint on w yielding w = — O^lg^gg with the flat 
prior on 6 2 . Results for the CMBall+B03 data set and the 
CMBall+B03+LSS+SNIa data set are given in Table 5. 

Figure 10 illustrates the degeneracy in the £l m — w plane 
that cannot be broken by CMB data alone and is only 
weakly broken with the addition of the LSS data (flat prior 
on b 2 g ). Application of a more restrictive Gaussian prior to 
b g for the LSS data or addition of the SNIa data breaks 
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Fig. 9. — Marginalized one-dimensional distributions for the 
dark matter equation of state parameter w. Weak priors im- 
posed on the base parameters are those outlined in Table 2. 
We also impose the prior —4 < w < 0. The t denotes the ap- 
plication of a Gaussian prior to b 2 g (with b g = 1± 10%). The 
nominal flat uniform prior on b g yields a slightly higher median 
value for w, driven by higher values of b g . Adding the SNIa 
data however, reduced the median value to —0.94 ± 0.1. 
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Fig. 10. — Constraints on w versus fl m for a flat ACDM model 
that allows the dark energy equation of state parameter, w, to differ 
from -1. Inner contours represent 68% likelihood regions and outer 
contours 95% likelihood regions. A more stringent Gaussian b 2 g prior 
(with b g = 1 ± 10%) or the addition of SNIa data is required to break 
the strong geometric degeneracy. 



It is a straightforward modification of CosmoMC to in- 
clude an additional isocurvature component in the MCMC 
chains, which uses CAMB to compute the isocurvature 
power spectra. Instead of using nj so as a basic parame- 
ter, we use two amplitude ratios for the two parameters 
that characterize our CDM isocurvature mode, following 
a suggestion of Kurki-Suonio et al. (2005): 

R 2 ee V iso (k 2 )/V s (k 2 ) , k 2 = 0.05M P c-\ (5) 

Ri = *W*i)/7M*i), h = 0.005MPC- 1 . 

Here V s (k) is the power in the primordial curvature per- 
turbation and Viso(k) is the power in the primordial CDM- 
photon entropy perturbation. The k 2 scale corresponds to 
I ~ 700 and k\ to I ~ 70. We adopt a uniform prior 
probability over the range to 20 for Ri, and over to 
100 for R 2 . The isocurvature spectral index ni so , defined 
by Viao(k) oc k nia ° , is now a derived parameter, expressible 
in terms of R\, R 2 and the adiabatic spectral index n s , 
defined by P s (k) cx k^' 1 , 

n iso =n s -l + \n(R 2 /R 1 )/ ln(k 2 /h) . (6) 

Two interesting limiting cases are the scale invariant nj SO = 
spectrum with high spatial correlation and the rij SO = 3 



"isocurvature seed" white noise spectrum with no spatial 
correlation 25 . 

The limits shown in Table 6 demonstrate that the large 
scale Ri, dominated by the WMAP data, is much better 
constrained at < 0.3 than the small scale R 2 < 2.3, which 
B03 probes. This translates into a preference for steeper 
nj so than the scale invariant value. Since for neither is 
there an indication of a non-zero value, just upper limits, 
the results are sensitive to the prior probabilities we assign 
them. Our choice of uniform prior for R\ and R 2 is conser- 
vative in that the upper limits decrease with other choices, 
e.g., one uniform in ln(i?^) (a non-informative prior), or 
one uniform in rii so and R 2 . The conservative choice actu- 
ally downgrades the probability of steep nj so . (The B03pol 
data by itself only limits R\ < 17 and R 2 < 22; the full 
B03 data gives R x < 1.8 and R 2 < 5.3.) 

Further constraints on isocurvature modes arise from 
LSS since the shape of the isocurvature matter power spec- 
trum differs in significant ways from the adiabatic one, but 
we do not consider those here. 

The strongest constraints come from the low I part 
of the spectrum 26 . However, spectra that are sig- 
nificantly steeper than inflation-motivated nearly-scale- 



as cosmic strings created in early universe phase transitions, which would contribute to a sub-dominant mass-energy content in the present 
epoch. Cosmic-defect-induced perturbations are greatly modified from their largely-uncorrelated initial state through gravity wave emission 
and have distinctively non-Gaussian features. By themselves they cannot explain the CMB data since the peaks and troughs are difficult to 
mimic, but they could be a sub-dominant component that future CMB data should be able to constrain. 

25 Such a spectrum is so steep that it must be regulated by a cutoff at high k >> k^. Physically this is typically the scale of the horizon when 
they are generated, but it is constrained by small scale structure information such as the allowed epoch of first star formation, which cannot 
be too early or else r could be far too large. We do not add this cutoff to our study since the CMB has a larger natural damping scale. A 
traditional seed case that has been considered is primordial black hole production. 

26 This is largely because of the isocurvature effect (Efstathiou & Bond 1986): to have no overall energy density perturbation on large scales, 
the entropy perturbation is carried almost entirely by the CDM or baryons, but when the equation of state changes from radiation to matter 
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invariant ones are still allowed by the data. To focus 
attention on the role played by the new, high I B03 re- 
sults, we now fix n; so at 3, the white noise 'seed' spec- 
trum, the limiting case in which the isocurvature pertur- 
bations when created were uncorrelated spatially. The 
large angular scales are highly suppressed and the isocur- 
vature peaks and troughs emerge looking somewhat like 
an ^-shifted version of the adiabatic spectrum. The two 
spectra then test at what level interleaved isocurvature 
peaks are allowed by the CMB data. Results are shown 
in Table 6. To relate the i?2 < 3.0 limit to a more in- 
tuitive expression of what the CMBall+B03+HST data 
set data allows, we note that over a bandpower in I from 

75 to 1400, C^ T(iso) /Cl T(s) ~ (0.005)i? 2 , hence the up- 
per limit corresponds to an allowed CMB contamination 
of this sub-dominant component of TT of only a few per- 
cent. Over a bandpower in £ from 150 to 1000 we find 

Cf (iso) /Cg E{s) ~ (0.008)i? 2 , for the allowed EE isocur- 
vature bandpower contamination. B03pol gives R2 < 58. 
The full B03 data set including TT gives R 2 < 9.5 27 . 

Since the isocurvature models adopted for these tests 
are not especially well-motivated physically, we have also 
chosen not to apply the LSS prior to our results. 

The n; so = 3,2 illustration allows us to conclude that 
even with the errors on the EE and TE data, there is ev- 
idence against the isocurvature shifted pattern over the 
adiabatic pattern and only restricted room for an inter- 
leaved peak pattern, at a level below 50%. This test dif- 
fers from the adiabatic-only peak/trough pattern shift us- 
ing B03pol Fig. 3 since there are no interleaved peaks and 
troughs in that case. Examination of the CAMB models ob- 
tained from the marginalized constraints in Table 6 reveals 
that the parameters chosen by CosmoMC adjust to make 
the adiabatic Cf pattern compensate for the isocurvature 
Cf° contamination. 

6. CONCLUSIONS 

The B03 data set does well at constraining the cosmolog- 
ical parameters of the standard ACDM model. The results 
are in good agreement with those derived from other CMB 
experiments, as is evident in Table 4. The parameter con- 
straints derived from the B03 data set in combination with 
the WMAP data are highly competitive with those from 
the CMBall data set. 

We have applied the Slosar-Seljak modification to the 
WMAP data which has the general effect of broadening 
slightly the WMAP ID likelihood curves with the largest 
impacts on the (WMAP alone) median values of r (<~ 0.3<r 
increase) and Cl m (~ 0.6a decrease) for the baseline model 
and n run (~ 0.5ct increase) for the baselinc+running in- 
dex model. Our graphical representation of the standard 
model parameters in Figure 2 illustrates the impact of the 
addition of the LSS data, which shifts the median values 
very little. Figure 2 shows the best estimates for the pa- 
rameters of the ACDM model from current CMB and LSS 
rcdshift survey data. 



Our analysis of five extensions to the standard model 
is summarized in Table 5. Intriguingly we found two 
cases, the running index and neutrino fraction, where 
the addition of data from higher multipole CMB exper- 
iments drove median values away from the conventional 
value. The evidence for a running index however, is 
slight (< 2a). Our neutrino mass limit (assuming three 
species of nearly degenerate mass) from CMB data alone 
is m v < 1.0 eV and from the combined CMB and LSS data 
set is m v < 0.40 eV, with no b g prior, and m v < 0.16 eV, 
with b g = 1.0 ±0.10. 

We have explored the sensitivity of our CM- 
Ball+B03+LSS results to the prior imposed on the galaxy 
bias factor b g and have found that only the neutrino mass 
results (above) and dark energy equation of state results 
are significantly impacted. For the dark energy equation 
of state we find that applying a more restrictive Gaussian 
prior to b 2 g gives w = -0.97 ± 0.12 with b g = 1.0 ± 0.10. 
Without any b g constraint we find w = — 0.65±0.16. How- 
ever, addition of the SNIa data, with the flat prior on b g , 

yields w = —0.94^°'^. The results are consistent with 
the dark energy being the cosmological constant. 

While the polarization data is not yet at the level of 
accuracy of the intensity data, cross checks of best fit pa- 
rameters from the B03pol data and B03TT data indicate 
consistent results. The consistency of the shape parameter 
9 determined from B03pol and from B03TT demonstrates 
that the peak and trough positions forecast by the spectra 
are in robust agreement. Isocurvature modes are begin- 
ning to be constrained by the current CMB polarization 
data and our upper limits and phcnomcnological discus- 
sion represent a good starting point for future analysis of 
these more complex models. The CMB polarization data 
is emerging but is not yet driving parameter determina- 
tion. We look forward to future higher precision data in 
which polarization data will play a larger role. 

We gratefully acknowledge support from the CIAR, 
CSA, and NSERC in Canada, ASI, University La Sapienza 
and PNRA in Italy, PPARC and the Leverhulme Trust 
in the UK, and NASA (awards NAG5-9251 and NAG5- 
12723) and NSF (awards OPP-9980654and OPP-0407592) 
in the USA. Additional support for detector development 
was provided by CIT and JPL. CBN acknowledges sup- 
port from a Sloan Foundation Fellowship, WCJ and TEM 
were partially supported by NASA GSRP Fellowships. 
Field, logistical, and flight support were supplied by USAP 
and NSBF; data recovery was particularly appreciated. 
This research used resources at NERSC, supported by 
the DOE under Contract No. DE-AC03-76SF00098, and 
the MacKenzie cluster at CITA, funded by the Canada 
Foundation for Innovation. We also thank the CASPUR 
(Romc-ITALY) computational facilities and the Applied 
Cluster Computing Technologies Group at the Jet Propul- 
sion Laboratory for computing time and technical support. 
Some of the results in this paper have been derived using 



dominance the photons carry the perturbation. The result is a rather dramatic amplification in the power over the Sachs- Wolfe power familiar 
in the adiabatic case by 36. This is why the nearly scale invariant case for isocurvature alone has been ruled out for so long from CMB data, 
so there was only room for it appearing at a sub-dominant level. 

27 The fi[ so — 2 case mimics even more the peak/trough patterns in except for the shift, so we tested that case as well. CMBall-l-B03+HST 
gives R2 < 2.7, B03 alone, but with TT, gives R2 < 6.8 and B03pol gives R2 < 41. Translation to the allowed contamination is done with the 

bandpower ratios Cg T(iso) /Cg T(s) ~ (0.007)^2, C% B(iso) /cf E(s) ~ (0.009)i? 2 . 
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Table 1 

B03 bandpowers. The lowest bandpowers of the TT spectrum (I < 375) are excluded when combining the B03 data with the 
WMAP results since the two spectra are signal dominated and therefore correlated. 



B03 Spectrum 


Multipole Range 


Number of Bands 


Reference 


TT 
TE 
EE & BB 


75 (375) <£< 1400 
150 < £ < 950 
150 < £ < 1000 


24 (18) 
9 
7 


Jones et al. (2005b) 
Piacentini et al. (2005) 
Montroy et al. (2005) 



Table 2 

List of weak priors imposed on baseline parameter set. Priors are uniform in the variable shown. 



Parameter 


Limits 


Parameter 


Limits 


n b h 2 
n c h 2 
e 

T 


0.005 - 0.1 
0.01 - 0.99 
0.5 - 10.0 
0.01 - 0.8 


n s 

ln[10 10 A s ] 
Age(Gyr) 
Ho 


0.5 - 1.5 
2.7-4.0 
10 - 20 
45 - 90 



Table 3 
The CMBall data set. 



Experiment 


Multipole Range 


Reference 


WMAP TT 
WMAP TE 
DASI TT 
VSA TT 
ACBAR TT 
MAXIMA TT 
CBI TT 


2 < £ < 899 
2 < £ < 512 
380 < £ < 800 
400 < £ < 1400 
400 < £ < 1950 
450 < £ < 1150 
750 < £ < 1670 


Hinshaw et al. (2003) 

Kogut et al. (2003) 
Halverson et al. (2002) 
Dickinson et al. (2004) 

Kuo et al. (2004) 
Hanany et al. (2000) 
Readhead et al. (2004a) 
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Table 4 

Marginalized parameter constraints for the baseline, 6 parameter, ACDM model. Parameter uncertainties represent the 68% 
confidence interval obtained by integrating the marginalized distributions. 95% confidence limits are quoted for the case of upper 
bounds. The following flat weak priors are imposed (as outlined in Table 2): 0.5 < n s < 1.5; 2.7 < ln(10 10 A s ) < 4.0; 
0.005 < Q, b h 2 < 0.1; 0.01 < Q, c h 2 < 0.99; 0.5 < 6 < 10.0; and 0.01 < r < 0.8. Additional weak priors restrict the age of the 
universe to lOGyr < age < 20Gyr and the expansion rate to 45 < Ho < 90. The CMBall data set is as given in Table 3. The LSS 
data consists of the galaxy power spectra from the 2dFGRS and SDSS redshift surveys. Our baseline CMBall+B03+LSS result 
is fairly insensitive to Sb g and we have chosen for this case the less restrictive fiat, uniform prior in b 2 g . The constraints from the 
B03pol data are in good agreement with the B03TT data, although some parameters constraints for the B03pol case are prior 
driven, eg. n B ,A B and Ho- B03 does not constrain r, but upper limits are given. The constraints from the various data set are 

consistent. 





B03pol 


B03TT 


B03 


WMAP 


WMAP 


CMBall 


CMBall 










+B03 


+B03 


+B03+LSS 


S\h 2 


0187+ ' 0058 

U.U-LO I _ .0058 


O.0219+° °° 3 3i 


0218+ ' 0029 
u.uzio_ 0030 


0245+ 0021 
u.uz4o_ 0019 


0241+ 0017 


U.VZM_ Q012 


n n997+ 000s 

U ' UZZ ' -0.0008 




n 1 rc + 0.047 
u - lo0 -0 .045 


0.129±S-~S 




m+0.016 


109+ 013 

u - ■ Luy -0 .013 


n i n« + D -010 
°'1° 6 -0.010 




8 


i.080tg;S~ 


i.uoi_ 010 


1 n55 +0010 
i.uoo_ 010 


i n4s+ 007 
l.U48_ 007 


1 049+ 005 

A .u*±y — 0.005 


1 045 + ° 004 
1 - u *°-0.004 




T 


< 0.66 


< 0.50 


< 0.49 


23+ 013 
u -"_0 .11 


22+ ' 11 
u -^ z -o.io 


17Q +0.090 


106+ 047 

U.J.UO_ 04g 


n B 


91+ 59 


n «K+ 010 

O-SD-0.10 


0.86_ 010 


1 02+ 06 


1 oi +0 - 06 

- L - UJ --0.05 


°- 98 -0.03 


95+ 02 

-0.02 


ln[10 10 A B ] 


3 4+°' 6 


3 4+°' 6 


3 4+°- 6 


3 4+°' 2 


q ,+0.2 
— 0.2 


o 9+0.2 
°- z -0.2 


, ,+0.1 

a - 1 ~o.i 


Oa 


o.59tg;5i 


o on+°- 19 

U.DJ_ 22 


n r r+0.17 

0.65_ 19 


o.76ig;g 7 7 


77 +0.06 
u - ' ' -0.06 


77 +o.05 

u - ' ' -0.05 


o.7oir 3 


Age(Gyr) 


i3.0ti;S 


jo r + 0.6 


13.4+g-f 


13.3«- 4 4 


, , O + 0.3 
lc, -°-0.3 


jo -+0.2 
lO.O_ 3 


13.6lg; 2 


O m 


41+ - 25 


0-37+oTo 


35+ ' 19 


24 +0.07 


n 90 + O.O6 
u - ZJ -0.06 


23+ - 05 


n on+0.03 
"■"-0.03 


OS 


1 2+°- 3 

-0.3 


o.98tg;i? 


0.971°,-? 


94+ 12 


o.92l° ;° 9 9 


n OO+0.06 


85+ - 05 
u -°°-0.05 


Z r e 


„+16.7 
J0 - z -17.7 


99 o + H.l 




19-611;? 


I9.3l 5 5 ;? 


16.5l 5 5 ;t 


12 6+ 3 ' 9 


Ho 


70 n +20 '° 

,u ' u -25.0 


69.2l 24 ;?, 


70.4±2 9 5-.4 


76.3± 4 3 2 7 


76.7l 3 3 9 3 


75.8l|? 


69.61 2 ; 4 
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Table 5 

Marginalized parameter constraints for five modifications of the baseline model. Parameter uncertainties represent the 68% 
confidence interval. For the case of upper or lower bounds 95% confidence limits are quoted. The following flat weak priors are 
imposed on the base 6 parameters (as outlined in Table 2): 0.5 < n a < 1.5; 2.7 < ln(10 1( M s ) < 4.0; 0.005 < Q. b h 2 < 0.1; 
0.01 < n c h 2 < 0.99; 0.5 < 9 < 10.0; and 0.01 < r < 0.8. Additional weak priors restrict the age of the universe to 
lOGyr < age < 20Gyr and the expansion rate to 45 < Ho < 90. We add, in turn, one parameter to the base set and impose the 
following prior on each: running index, —0.3 < n run < 0.3; curvature, —0.3 < Q.k < 0.3 ; amplitude ratio, < A t /A s < 20; 
neutrino fraction of dark matter, < /„ < 1; and dark energy equation of state, — 4 < w < 0. The CMBall data set is as given in 
Table 3. The LSS data consists of the matter power spectra from the 2dFGRS and SDSS redshift surveys. We have marginalized 
the galaxy bias factor b g assuming a uniform flat prior in b g . For the Base + w case only we add in both the LSS and SNIa data 
which gives a better constrained result than that which is obtained from adding LSS data alone. We include the baseline model 

results for comparison. 





Base-f- Running Index 


B as e-f- Curvature 




CMBall 
+B03 


CMBall 
+LSS+B03 


CMBall 
+B03 


CMBall 
+LSS+B03 


n b h 2 

6 

T 

n s 

ln[10 10 A s ] 

Oa 
Ago(Gyr) 

^8 

Z re 

Ho 
n run 


104+ 015 

u ' lu *-0 .013 

1 n4s+ 005 

UJ1 _0.14 

95+ 06 

q c+0.2 
J - o -0.3 

i3.3l° ; 4 4 

0.22±°$ 
24.7+_U 

7S ,+11.9 

78.1_ 8 7 

-n n72+ - 036 

U - U,Z -0.036 


0219+ 0010 
u.uzi»_ ooio 

125+ 007 

1 044+ 004 

1,um -0 .004 

151+ 025 

U.1UJ._ o3 6 

90+ 04 

U.»U_ 04 
q 9 + 0.1 

a - z -o.i 

o.67t° :°„ 4 4 
i3.7i° : 2 2 

n o, + 0.04 

0.89t o:°06 

i6.8i 5 5 ;°„ 

n 051+ ' 027 

U.UOl_ 026 


0231+ 0010 

U.UZOl_ 0013 

106+ 011 

U.1UD_ on 

1 045+ 005 

-L.U<iO_ 005 
0.178_ 092 
97+ 03 

U.VI _ 04 
, „ + 0.2 
°- z -0.2 

67+ ' 12 

-0.13 

14 9+ 13 

n o 7 +0.17 
"■•>' -0.16 
n oi +0.06 
U - ai -0.06 

16.6+^ 7 

62.B±i*;i 

-0.037t° ;°o 33 9 


n n997+ -°oo 9 

u.uzz 1 -o.ooiO 

111+0.008 
u - lii -0.0O8 

1 044+oo° 4 

1,um -0 .004 

139+ 021 

u - 1Ja -0 .129 

96+ 02 

u ' MD -0.03 
n 1+0.1 
31 -0.1 

67+ 05 

-0.04 

14 8 +0 ' 7 

n oc + 0.05 
n o 9 + 0. 06 

i4.7i 4 4 ; 7 

6i. 7 ; 2 ; t 

-0.027t :°oi 6 6 




Base-j-Tensor Modes 


Base-j-Massive Neutrinos 




CMBall 
+B03 


CMBall 
+LSS+B03 


CMBall 
+B03 


CMBall 
+LSS+B03 


n b h 2 
n c h 2 

6 

T 

n B 

ln[10 10 A s ] 

Oa 
Age(Gyr) 

^8 

Zre 

Ho 

A^ 1 A 
fit / ^-s 

fv 


0246+ 0013 

U.UZ+D_ 0013 

0957+ 009 " 

i r>A£+ - 004 
1.048_ 004 

162+ 072 

U.lO^_ 070 

1 02+ - 04 

o ,+0.1 

^■'■-o.i 

n HO +0.03 

, , ,+0.3 

(1 iai+0.036 
U.181_ 034 

0.77tS;gJ 
15.2± 4 4 j 

> 73.2 

< 0.71 


n n9q9+ 0009 

U.U^O^_ QQog 
n , , 7 +0.006 
-0.006 

1 046+ 004 

l.U+D_ 004 

104+ 049 
u.iui_ 04g 

97+ 02 

u ' y ' -0.02 

,+0.1 

a - 1 -o.i 
0.72t :°o 3 3 

i3.5i° ; 2 2 

o.28i° :°„ 3 3 

0.84t° :°o 5 5 
12 2+ 4 ' 
71 4+ 2s 

1 ^'^-2.9 

< 0.36 


0226+ 0014 

U.UZZO_ .0014 

120+ 014 

U.1ZU_ 014 

1 047+0 005 

J-'"*' -0.005 

164+ 094 

U.1D1_ 088 

95+ 04 

U.»O_ 04 
o 9 + 0.2 
°- z -0.2 

64+ 11 

U ' D *-0.11 

i4.ii : 4 4 

o.eoiri 3 
16. 6 ; 6 5 ; 2 

64.9t 7 ; 4 
< 0.21 


0224+ - 0008 

u.uzz^ o.0009 

126+ 007 

U.1ZO_ 007 

1 045 + 00 4 
'■ u,J -0 .004 

n ins+o- 049 
U.1U8_ 047 

95+0- 02 

U ' MO -0.02 
n 1+0.1 

•'■'--o.i 
64+ 06 

u '°^-0.05 

13 9 +0 ' 2 
io.y_ 2 

3R+O O 5 
0^t°oZ 

i3.oi 4 4 ;„ 

64.8i 3 ; 9 
< 0.09 




Base-j-w 


Baseline 




CMBall 
+B03 


CMBall 
+LSS+B03+SNIa 


CMBall 
+B03 


CMBall 
+LSS+B03 


n b h 2 
n c h 2 

9 

T 

n B 

ln[10 10 A s ] 

Oa 
Ago(Gyr) 

n m 

z re 
Ho 
w 


0234+ 0013 

u.uzoi_ 0013 

106+ 011 

U.1UD_ on 

1 046+ 005 

l.U+D_ 005 

161+ 078 

u.lul_o.077 
, „ + 0.1 

°- z -o.i 
71+ - 10 

13 7+ - 4 

LO - ' -0.4 

29+ ' 11 
u.z»_ i 

n -t^ + 0.15 
°- 76 -0.15 

15 9+ 50 

-o.86±°-.iS 


0229+ 0009 

U.U^ZJ_ 0009 

n 117+0.007 

-0.008 
1 045 + 004 
±.U1O_ 004 

0-H8 + n 5 55 

96+ 02 

"• au -0.02 

o ,+0.1 

•J^-o.i 
70+0- 02 

u - ,u -0.02 

13 6 +0 ' 2 

lo.u_ 2 

30+0' 02 
U - 3U _0.02 

n OO+0.06 

U -8^_0.06 

13 3+ 4 ' 2 

lo.o_ 4 2 

68.6; 2 ; 

_0 94+0- 094 

u,y4: — 0.096 


n n9q-?+o o° 13 
u.uzoo_ 0012 

106+ 010 

U.1UD_ oio 

1 045+0- 004 

l.WiO_ () 004 

170+ 090 

u ±l "-0.081 

n qs+O- 04 
u - y 8-o,03 

, 9 + 0.2 
°-^-0.2 
77 +0.05 
u - ' ' -0.05 

23+ ' 05 

u - ZJ -0.05 

°.83±r„ 6 6 

16-5i 5 5 4 
75.8t 5 5 « 


n n997+ -000 s 

U.UZZ(_ 0008 

120+0- 005 

U.1ZU_ 005 

1 045 + 00 4 
'■ u,J -0 .004 

106+ 047 

U.1UU_ 048 

95+ 02 

u ' aJ -0.02 

3 1+ - 1 
^•^-o.i 

70+ 03 

u - ,u -0.03 

13 6 +0 ' 2 
io.u_ 2 

30+ 03 
u .°°_0.05 

i2. 6 ; 3 - 9 

69.61 24 
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Table 6 

Marginalized parameter constraints for a model which includes both (dominant) adiabatic and (sub-dominant) isocurvature 
modes. Parameter uncertainties represent the 68% confidence interval. Upper bounds are 95% confidence limits. The flat weak 
priors are imposed on the base 6 parameters are as outlined in Table 2. The CMBall data set is defined in Table 3. We include 

the baseline model result (with the more stringent HST prior) for comparison. Wc consider two parameterizations for the 
isocurvature model. For the first (column two) we add two parameters to our basic six: R 2 = 7 , i B o(fe) /V B {k,2) , with pivot scale 
k 2 = 0.05Mpc _1 ; and i?i = V iso (ki) /V a (ki) with fci = 0.005Mpc _1 . We impose the priors < Ri < 20 and < R 2 < 100. For 
this case the isocurvature spectral index, rii so , is a derived parameter. We also consider the "white isocurvature" case (column 
three) where we fix rii so = 3 and allow the amplitude ratio R2 = (Ai ao /A s ) to lie anywhere between and 100. 





Baseline 


Adiabatic + Iso 


Adiabatic + White Iso 




CMBall 
+B03+HST 


CMBall 
+B03+HST 


CMBall 
+B03+HST 


n b h 2 
n c h 2 

6 

T 

n a 
Ri 
Ri 

«i B 

ln[10 10 A s ] 

n A 

Ago(Gyr) 
H 


0230+ 0010 
u.uzou_ 0010 

U.1U8_ 008 

1 045+ 004 

1 - u,J -l) .004 
1C1 +0.069 
ulol -0.067 

97+ 03 

u.y < _Q.03 
,+0.1 

76+ ' 04 

lO.O_ 2 

o.24i ;° t 

74.4+. 4 ; , 


U-U^*' -0.0014 
„ ,„o + 0.009 
u - luo -0 .009 

i.uoi_ 005 
n ifiq+0 071 

1 01+ 004 

1 ' ul -0.04 

< 0.28 

< 2.3 

1 1+0-8 
q 9 + 0.1 

J - z -o.i 
n sn+ - 04 

13 l+°' 3 

0-20i o o ;° o 4 4 

i5.s; 4 4 - 4 4 

80.3t 5 5 6 5 


n239 + 0012 
u.umo_ .0011 

q 10 y+0.009 

1 046+ 005 

15fi + 071 

97+ 03 

u - y ' —0.03 

< 3.0 
3.0(/izeci) 

, + 0.1 

77 +0.04 
u - ' ' -0.04 

13 4 +0 ' 2 
0-23i ;° 4 4 

i5-6i 4 4 ; 8 6 

75.9t 4 4 ; 3 



